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ABSTRACT 

For growth of a vicinal face at constant velocity, the effect of anisotropic interface kinetics on mor- 
phological stability is calculated for a binary alloy. The dependence of the interface kinetic coefficient on 
crystallographic orientation is based on the motion and density of steps. Anisotropic kinetics give rise to 
traveling waves along the crystal-melt interface, and can lead to a significant enhancement of morphological 
stability. The stability enhancement increases as the orientation approaches a singular orientation and as 
the solidification velocity increases. Shear flows interact with the traveling waves and, depending on the 
direction of the flow, may either stabilize or destabilize the interface. Specific calculations are carried out 
for germanium-silicon alloys. 

INTRODUCTION 

During the directional solidification of a binary alloy solute inhomogeneities can arise from both fluid flow 
and morphological instability. In microgravity buoyancy-driven fluid flow is reduced, and experiments to 
study the evolution of morphological patterns without the interference of fluid flow may be possible. We are 
carrying out theoretical studies of the interaction of fluid flow with the crystal-melt interface [1-3]. Included 
in this research are: (1) calculations of cellular morphologies in the absence of fluid flow; (2) evaluation of the 
Seebeck voltage for cellular interfaces as a method for monitoring interface morphology in metallic alloys; 
(3) linear stability analyses of coupled interfacial and convective instabilities; (4) calculations of the effects 
of time-dependent gravitational accelerations (g-jitter) on fluid flow during directional solidification. This 
ground based research will focus on providing theoretical interpretation and guidance for a series of space 
experiments to be carried out by J. J. Favier, R. Abbaschian, and colleagues on tin-bismuth alloys using 
the MEPHISTO apparatus and by K. Leonartz and colleagues on succinonitrile-acetone alloys. Bismuth 
typically grows from the melt with facets, and in this paper we will describe recent studies of morphological 
stability when the interface attachment kinetics is highly anisotropic. 

During solidification of a binary alloy, the crystal-melt interface may become morphologically unstable 
[4, 5] leading to cellular or dendritic growth. Here, we consider the effect of anisotropy of interface kinetics 
on morphological instability during growth of a binary alloy at constant velocity. We treat growth in a 
direction near a singular orientation, and use a model of kinetic anisotropy based on step motion. The 
effect of anisotropy of surface tension and interface kinetics has previously been treated [6] in a quasi-static 
approximation to the diffusion field; kinetic anisotropy gave rise to traveling waves along the crystal-melt 
interface. In this article, we keep the full time-dependence of the diffusion field and show that the traveling 
waves resulting from kinetic anisotropy can provide a significant enhancement of morphological stability, 
especially for orientations near the singular orientation and for large growth velocities. We have recently 
obtained similar results for constant velocity growth from a supersaturated solution [7]. 

In solution growth, it is well known that a solution flowing above a vicinal face of a crystal has a strong 
influence on morphological stability [8-10]. If the flow is opposite to the direction of the step motion, there 
is a stabilization of the interface with respect to Btep bunching. In the absence of a shear flow, there is an 
apparent flow of liquid opposite to the step motion (in a reference frame attached to the moving steps) and 
this stabilizes the interface, preventing step bunching and the formation of macrosteps [7]. 

For melt growth, interface stability near a vicinal face was the subject of an earlier treatment by Yuferev 
[11] in which he concluded that interface stability is enhanced for very small deviations from the singular 
orientation. Also, Kolesnikova and Yuferev [12] have carried out weakly nonlinear calculations in the frozen 
temperature approximation, i.e., equal thermal properties and vanishing latent heat. A strong interaction 
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between a shear flow parallel to the crystal-melt interface and anisotropy of kinetics has recently been found 

[ 10 ]. 

In this article, we carry out a linear treatment of interface stability during constant velocity solidifica- 
tion of a binary alloy. We allow for an interface kinetic coefficient that depends on the crystallographic 
orientation of the crystal-melt interface. We carry out specific calculations for germanium-silicon alloys for 
orientations in the vicinity of a singular orientation, which occur in this alloy [13]. 


THEORY 


We consider directional solidification of a binary alloy at constant velocity Y in the 2 -direction and treat 
the morphological stability of an initially planar crystal-melt interface. We choose an (x, y, z) coordinate 
system (moving with the macroscopically planar interface) such that the crystal-melt interface is described 
by 2 = h(x,t), where t is the time. Further, we consider a two-dimensional problem and assume that all 
quantities are independent of the ^-coordinate. 

The temperature fields, Tjj(x,z,t) and Ts(x,z,t), in the melt and crystal and the solute field in the 
melt, Ci,(x,z,t), satisfy the partial differential equations 
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where kl and k$ are the thermal diffusivities in the melt and crystal, respectively, and D is the solute 
diffusivity; we neglect diffusion in the solid. We also assume that there is no fluid flow in the melt and that 
the thermophysical properties are independent of temperature and solute concentration. These transport 
equations and assumptions are most appropriate for dilute alloys. 

The boundary conditions at the crystal-melt interface are 


(V ■ n)(l - k)C L = -DVC l ■ n, (4) 

L v (V • n) = ( k s VT s - k L VT L ) • n, (5) 

T s = T l , ' (6) 

V • n = P(p)(T e — n), (7) 

where 

T e = Tm + mCx — TmT/C, (8) 


V is the solidification velocity vector, n is the unit normal to the interface, k is the segregation coefficient, L v 
is the latent heat per unit volume, ks and fen are the thermal conductivities of crystal and melt, respectively, 
f. ?(p ) iB the kinetic coefficient, Tm is the melting point, m is the liquidus slope, T = j/L v , 7 is the surface 
tension, and fC is the mean curvature of the crystal-melt interface. We will assume that (3(p ) = P,t\p\, with 
p = tanfl, where 9 measures the deviation of the interface from a singular orientation. We will denote the 
orientation of the planar interface relative to the singular orientation by p, and consider that p > O.We fix the 
far-field boundary conditions by specifying the bulk alloy concentration C x and the temperature gradient 
in the melt at the crystal-melt interface. We will perform stability calculations for germanium-silicon alloys; 
the required properties are given in reference [14]. 

In [7, 10, 14] the morphological stability of a solid-liquid vicinal interface is treated using linearized theory 
in two dimensions. The perturbation to the planar interface (2 = 0) with crystallographic orientation p has 
the form 

2 oc exp(crt + ikt.x). (9) 
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Here k x is the wavenumber and cr is the complex temporal factor, 

a = < 7 , + itTi. (10) 

The interface is unstable if oy > 0. The dispersion relation for a as a function of k x can be obtained 
analytically. In the thermal steady-state approximation, i.e., — » oo and <65 — » 00, it reduces to 

ar/V = {-<? - T M Tkl + ik x V/(fi lt p 2 ) + mG e [a - V/D]/[a - (1 - k)V/D}} 

{L v V/{2 kk x ) + (V/p tt p) + mG e /[a - (1 - k)V/D]}~\ (11) 

where 

a = (V-/2 D) + y/(V/2 D)> + k* + (v/D), (12) 

<5 = ( ksGs + ki J Gi l )/{2k) } and k = (ks -I- k£)j2. In a previous treatment [6] of the effect of anisotropic 
kinetics on morphological stability, a quasi-static approximation was made for the diffusion field which is 
equivalent to setting <r = 0 in the definition of a. In this limit, the above expression agrees with eq. (13) of 
reference [6]. The dispersion relation can be written as a cubic polynomial with complex coefficients in the 
variable a — (V /2 D)\ this has the advantage that all roots of the nonlinear equation can be determined. In 
the thermal steady-state approximation, we have numerically solved the cubic polynomial to determine the 
stability conditions. We have also numerically solved the ordinary differential equations as a boundary value 
problem to determine the stability conditions; this does not require the thermal steady-state approximation. 

RESULTS 

The dispersion relationship ay = ay (y, G, k x ) includes, as independent parameters, the normal 

growth rate V, the temperature gradient G, the impurity concentration and the slope p; stability requires 
that ay < 0 for all values of k x . In principle, the slope, p, may be fixed independently of the other processing 
variables by cutting the crystal to prepare an initial interface; however, only very small deviations (less than 
a tenth of a degree) from a singular interface make sense because large deviations are indistinguishable from 
a rough interface. Even such a prepared interface may be replaced (wedged out) by the nearest singular face 
in the course of sufficiently long growth. Under the typical conditions of faceted growth, the vicinal slope 
p on a facet is, like V , dependent on supercooling and is determined either by a dislocation step source or 
by two-dimensional nucleation on the most supercooled region of the facet. In these cases, the slope p is 
a function of the growth rate V which is determined by the conditions of crystal pulling (Czochralski) or 
crucible movement through a temperature gradient (Bridgman and related techniques). In the following, 
we consider stability conditions for both an independently given slope and for a slope determined by V, i.e., 
p = p(V), for dislocation assisted growth and for growth by two-dimensional nucleation. 

Figure 1 presents the real part of the temporal factor oy i.e., oy, as a function of the wavenumber k x in 
dimensionless form for the processing conditions V = 10 pm/s, Gi = 50 K/cm, and C*, = 5 at%. The 
upper plot with p = 0.002 is for isotropic kinetics. The lower two plots are for anisotropic kinetics with p 
= 0.002 and 0.0015. Stabilization due to anisotropy is evident. Also, for p < 0.0015 the interface becomes 
absolutely stable, i.e., stable against perturbations for all k x . Stability at Dk x /V < 5 is common to all the 
curves; in this region of k x the perturbation wavelength 2ir/k x > D/V, i.e., the diffusion boundary layer 
follows the perturbed interface and thus there is less driving force for instability. The steep drop in oy at 
high wavenumbers is associated with capillary stabilization, while the kinetic stabilization determines the 
overall level of a T at lower k x and thus the height of the maximum. 

We suppose that the unperturbed slope p is determined by a screw dislocation with the elementary 
Burger’s vector intersecting the face under consideration. We also assume that the supercooling, AT, and 
thus the value of p, along the face is constant with AT given by V /{/3,tp). Letting j, t be the linear energy 
of a step, the growth rate of the face is given by [15] 

= AT» (13) 

197..TM 

and the slope is 

p=[h 3 L v 9/(19 7 , t T M 0, t )]5 (14) 
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where h is the step height. It was shown by Voronkov [15] that the step fluctuations at the melting point 
decrease the step free energy for (111) Si and Ge by a factor of three or four; therefore, for the linear stability 
calculations for Ge, we assume that y 4 * = 7 / 1 / 4 . 

If the crystal face is free of screw dislocations or the supercooling is so high that step generation by 
two-dimensional nucleation is more effective than by the dislocation mechanism, the growth rate is [16] 

V = h(v 2 J)5 (15) 

where the nucleation rate J (cm -2 s) is obtained from 

J=[hk B T$ I AT/(w 2 L v )]%ql3, t exp[-Trhy 2 /(k B ATL v )]. (16) 

Here q is the density of surface sites on which nucleation is possible, u is the atomic volume, and k B is 
Boltzmann’s constant. Thus, the growth rate is 

V = {h 7 q 2 /u>)\p. t AT *[k B Tjl,/(uL v )]Z exp[-xhy 2 /(3k B ATL v )}. (17) 

and the slope p is given by A 

p— h(J/v)3 = (h 7 q 2 /u)6 ^ tB 'T /^)T e xp[—xhy 2 /(3k B ATL v )]. (18) 

ATe 

Thus eq. (17) and eq. (18) give the dependence p(V) in parametric form with AT the parameter. In 
evaluating these expressions, we make the approximation that h 7 q 2 /w = 1 . 

In Fig. 2, we present a stability diagram in terms of the critical solute concentration as a function of 
the ratio of the liquid temperature gradient and the growth velocity for G B = 50 K/cm. The interface is 
stable to small perturbations for concentrations lying below the given curves. The calculations have been 
carried out by numerically solving the differential equations. The solid curve is the Mullins and Sekerka 
[4] result which for isotropic kinetics does not depend on the kinetic coefficient [ 6 ]. The dashed and dotted 
curves are calculated from the present model of anisotropic kinetics with p = 0.01 and 0 . 001 , respectively. 
Anisotropic kinetics causes a substantial enhancement of the region of stability. 

The remaining two curves in Fig. 2 are for dislocation and two-dimensional nucleation controlled growth 
for which p is a function of velocity. For the dislocation controlled growth, V = 8.7 X 10 4 p 2 cm/s. In the 
case of two-dimensional nucleation, eq. (17) and eq. (18) have been evaluated numerically to obtain V as a 
function of p. 

We note that as p becomes Bmall, the concentration at the onset of instability becomes very large where 
other effects, which we have not taken into account, such as fluid flow and the concentration dependence of 
the distribution coefficient, liquidus slope, and diffusion coefficient, may be important. However, it is clear 
from the results that as a singular orientation is approached, anisotropic kinetics has a large stabilizing 
effect. 

As p becomes very small or V becomes very large, the thermal steady-state approximation is no longer 
valid. For example, in the thermal steady-state approximation, for p = 0.001 and V = 0.001 and 0.01 cm/s, 
the critical concentrations are 7.8847 and 5.4579 at%, respectively; whereas the numerical values are 7.8852 
and 7.6689 at%. Thus for p = 0.001, the thermal steady-state approximation underestimates the critical 
concentration. The thermal steady-state approximation fails because the wavespeed is sufficiently large that 
k 2 is no longer much greater than o»/ks. 

Near a singular orientation, due to anisotropy the value of TmT may be much less than the value 
used in our calculations [14]; therefore, we have also carried out a few calculations with TmT reduced 
by a factor of 100. For p = 0.001 and V ranging from 0.001 to 0.1 cm/s, the critical concentrations are 
essentially unchanged. Similarly, for dislocation controlled growth, the results for the critical concentrations 
are again unchanged. Thus, in the strongly kinetically controlled regime the value of the capillary constant 
is unimportant. Capillarity is still important at large wavenumbers, but in the kinetically controlled regime 
the onset of instability occurs at relatively small wavenumbers. 
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CONCLUSIONS 


In the calculations described above, we have assumed that the melt is stagnant. It is well known that 
fluid flow can alter the conditions for stability [10, 17]. For an ideal fluid and anisotropic kinetics, Chernov 

[10] has shown that a flow parallel to the interface and opposite in direction to the step motion stabilizes 
the interface while a flow in the direction of the step motion destabilizes the interface. We are carrying 
out numerical calculations of the linear stability of a crystal-melt interface in the presence of a parallel 
shear flow. The method is similar to that previously used for isotropic kinetics [18]; in addition, we have 
developed a pseudospectral algorithm to treat the problem. Qualitatively, the results are in agreement with 
Chernov [10] for an ideal fluid. 
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Figure 1: The real part of the temporal growth rate as a function of the spatial wavenumber for three values 
of the crystallographic orientation with respect to the growth direction for p = 0.002 and 0.0015. The upper 
plot is for isotropic kinetics, while the two lower plots are for anisotropic kinetics. 



Figure 2: The silicon concentration at the onset of instability during directional solidification of germanium- 
silicon alloys as a function of the liquid temperature gradient divided by the growth rate with G/> = 50 
K/cm. The interface is morphologically stable for concentration values below the curves. The solid curve 
is based on the Mullins and Sekerka theory with isotropic kinetics; the remaining curves are for anisotropic 
kinetics. The dashed and dotted curves correspond to p — 0.01 and 0.001, respectively. The curves labeled 
“nucl” and “disl” correspond to nucleation and dislocation controlled growth for which p is a function of V. 
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